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ABSTRACT 

We present the detection and modeling of more than 70 far-IR pure rotational lines of water vapor, 
including the 18 and 17 isotopologues, towards Orion KL. Observations were performed with the 
Long Wavelength Spectrometer Fabry-Perot (LWS/FP; A/AA ^6800-9700) on board the Infrared 
Space Observatory (ISO) between ~43 and 197 /im. The water line profiles evolve from P-Cygni type 
profiles (even for the Hj O lines) to pure emission at wavelengths above ~100 /im. We find that most 
of the water emission/absorption arises from an extended flow of gas expanding at 25±5 km s~ 4 . 
Non-local radiative transfer models show that much of the water excitation and line profile formation 
is driven by the dust continuum emission. The derived beam averaged water abundance is 2-3 xlCP 5 . 
The inferred gas temperature Tfc=80-100 K suggests that: (i) water could have been formed in 
the plateau by gas phase neutral-neutral reactions with activation barriers if the gas was previously 
heated (e.g. by shocks) to >500 K and/or (ii) H2O formation in the outflow is dominated by in-situ 
evaporation of grain water-ice mantles and/or (iii) H2O was formed in the innermost and warmer 
regions (e.g. the hot core) and was swept up in «1000 yr, the dynamical timescale of the outflow. 
Subject headings: infrared: ISM: lines and bands — ISM: individual (Orion) — ISM: molecules — line: 
identification — radiative transfer 



1. INTRODUCTION 

Star forming regions are associated with violent phe- 
nomena such as cloud collapse, molecular outflows and 
related shocked regions. Under these conditions the neu- 
tral g as is warm and water vapo r is predicted to be abun- 
dant l)Draine fc Robergelll982D and to play a dominan t 
role in the thermal balance jNeufeld fe Kaufmanlll993|) . 
Unfortunately, ground-based observations of H2O are 
difficult and so the determination of water column densi- 
ties is not straightforward. Nevertheless, ground-based 
observations of water maser lines have been performed 
toward Orion KL. F rom VLBI observat ions of the 6i6~ 
523 line at ~22 GHz. lGenzel et aTll)1981|) determined the 
kinematics and the expansion velocity for the so-called 
low velocity outflow (~18±2 km s _1 ). The widespread 
nature of water vapor has been probed with maps of th e 
3i 3 -2 20 line at -183 GHz l|Cernicharo et al.lll990lll99^ . 
This was the first time that the water abundance was 
estimated in the different large scale components around 
Orion KL. In addition, the high excitation conditions of 
the plateau gas (a mixture of outflows, shocks and inter- 
actions with the ambient cloud) has been r evealed by ob- 
serva tions of the 5i 5~4 2 9 line at ~325 GHz ijMenten et alJ 
1990: iCernicharo et al.lll999b|) . Due to their maser na- 
ture, the water abundance determination from these lines 
is quite involved. Even the observation of the H^O 3i3~ 
2 2 o line at ^203 GHz gives only poor estimates due to 
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the overlap with other molecular lines l)Jaca et al.l ll988). 
The HDO species has been also detected and modele d 
towards Orion KL ij.Iaca et al J 1 199(1 iPardo et alJl2001^ . 
but the specific gas and dust chemistry has to be taken 
into ac count to derive y(H 2 Q). 

ISO (Ke ssler et al.lll996|) spectrometers have provided 
the opportunity to observe many IR H2O lines towards 
bright sources such as Orion KL. In particu lar, Harwitt 
et al (1998) presented eight ISO/LWS/FP ijClegg et a.U 
1996) water lines between ~71 and ~125 ^m involving 
energy levels between 300 and 800 K. They estimated 
x(H 2 0)~5x 10 -4 ; however, the interpretation of these 
lines should include radiative pumping from IR dust pho- 
tons. A larger set of weaker IR lines, including those of 
H^O and H^O, is needed to minimize the opacity ef- 
fects always associated with HJ, 6 lines. Orion KL has 
also been targeted with the Short Wavelength Spectrom- 
eter (SWS) tracing a smaller region than that probed b y 
the LWS beam below ~45 ^m ijvan Dishoeck et alT l998i. 
Nineteen water absorption li nes with energies between 
200 and 750 K were detected l)Wright et alJl2000j) . These 
authors suggested that mid-IR water lines arise from the 
low-velocity outflow and estimated x(H20)=2-5x 10 -4 
(for an assumed gas temperature of 200-350 K). Most re- 
cently, the lio-loi lines of H^O and H 48 O have been ob- 
served with SWAS (Y(Hj, 6 Q)-3.5xlQ- 4 ; iMelnick et all 
(2000) :lSnell et alJ lEfflOTO') and ODIN (%(H| 6 O)~10" 5 - 
10~ 4 ; Olofsso n et a 1J (12003^ . In these observations, the 
H^O ground-state line shows widespread emission over 
a ~5' region. 

In this letter we present all water lines detected in the 
first far-IR line survey towards Orion KL carried out 
with the ISO/LWS spectrometer (Lerate et al. 2006), 
and the radiative transfer models that fit the water data 
set (involving levels up to 2000 K). 
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Fig. 1. — Summary of the 70 far-IR water lines detected by the ISO/LWS-FP towards Orion KL. The ordinate corresponds to the 
continuum normalized flux and the abscissa to the velocity (in km s~ 1 ). Transition rotational numbers, rest wavelengths (in fim), upper 
level energies (in K) and intrinsic line strengths are shown in each box. The central inset shows the ISO/LWS-grating observations towards 
Orion at a resolution of A/AA ~ 300. The ordinate corresponds to the absolute flux (in W cm -2 (im _1 ) and the abscissa to the wavelength 
(in (im). Main molecular features between ~160 and ~197 (im arc labelled in a zoom to the grating spectrum. The full radiative transfer 
model for the far-IR continuum and water line spectrum is also shown in the upper inset. 
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2. OBSERVATIONS AND DATA REDUCTION 

All water lines within the range of the ISO/LWS (-43- 
197/im) were observed using long integrations with the 
FP spectrometer (L04 AOT) which provides the largest 
spectral resolution for the LWS instrument. Preliminary 
results were presented bv iCernicharo et all l|1999a|L In 
addition, a complete less sensitive (L03 AOT) far-IR line 
survey of Orion has been carried out. Adding both data 
sets more than 70 water lines have been detected (Lerate 
et al. 2006). The LWS circular aperture size is ~80", 
although it slightly depends on the particular LWS de- 
tector. In its FP mode the instrumental response is close 
to a broad-wing Lorentzian with a spectral resolution of 
A/AA -^6800-9700. Processing of the water lines from 
AOT L03 was carried out using the Offline Processing 
(OLP) pipeline and the LWS Interactive Analysis (LIA) 
package version 10. AOT L04 spectra were analyzed us- 
ing the ISO Spectral Analysis Package (ISAP 6 ). The full 
description of the complex data calibration and reduction 
process and associated target dedicated time numbers, 
and all the tabulated line observations and spectroscopic 
parameters are given by Lerate et al. (2006). 

3. RESULTS AND DISCUSSION 

A summary of the resulting water lines is shown in 
Fig-E] From these far-IR observations it is clear that wa- 
ter lines show a complicated behavior. For wavelengths 
above -TOO /im, H3, 6 lines are observed in emission. 
However, for shorter wavelengths, lines arising from en- 
ergy levels below -4000 K and with large line-strengths, 
show P-Cygni profiles with emission covering a large ve- 
locity range. However, those with weak line-strengths or 
arising from higher energy levels, are observed in pure 
emission. H^O lines also show P-Cygni type profiles 
below -400 ^m. In the H3, 8 case, the absorption com- 
ponent is deeper than in the analogous H^O transition. 
Since optical depth effects are much less severe in H^O 
lines, the associated P-Cygni type profiles trace the main 
origin of far-IR water lines toward Orion KL, i.e. an ex- 
tended outflow. Pure emission H2O lines peak around 
v ~20 km s _1 (the source Vlsr is ~9 km s -1 ; Scoville et 
al. 1993) but it is likely that the most opaque lines are 
redshifted due to self-absorption. On the other hand, 
water lines detected below ~50 /im are observed in pure 
absorption with a velocity peak of v ~ —10 km s _1 (the 
same applies t o mos t water lines observed by ISO/SWS; 
IWright et all JSoOflh l. 

The turnover point between absorption and emission 
lines is an important clue to interpret this large data set 
and to determine the relations between continuum+line 
opacity, line strengths, spatial distribution of gas, and 
physical conditions. ISO observations clearly show that 
most of the water vapor detected in the IR arises from 
a flow of gas expanding at 25±5 km s — 1 . The in- 
ferred expansion velocity is consistent with the low- 
velocity outflow origin ally revealed by -^22 GHz H2O 
maser motions ijGenzel et al.lll98lj) . but a contribution 
from the extended high- velocity outflow could be present 
ijCernicharo et alH 9941. However, no high velocity line- 
wing emission is observed at ISO's sensitivity and S/N. 

6 ISAP is a joint development by the LWS and SWS Instruments 
Teams and Data Centers. Constributing institutes are CESR, IAS, 
IPAC, MPE, RAL and SRON. 



Similar conclusi ons have been found for the OH excited 
rotational lines ijGoicoechea et a l. 2006). 

The main problem with modeling ISO data is the lim- 
ited spatial resolution, which makes difficult to constrain 
the size and origin of the water region. Besides, any 
detailed fit to the data requires a detailed knowledge 
of the geometry and of the relative distribution of dust 
continuum and water lines. Fortunately, high angular 
resolution maps of water at -483 and ^325 GHz lines 
do reveal the spatial and velocity distri bution of water 
in Orion ijCernicharo et al1ll994l Il99 9b1. In particular, 
Cernicharo and coworkers detected at least four differ- 
ent water components: the extended ridge (extended 
quiescent gas), the plateau (including the high- and 
low-velocity outflows) and the very narrow and strong 
features associ ated with the small water bullets observed 
at ~22 GHz ijGenzel et all 11981(1 . Besides, the newly 
detected far-IR water lines associated with the highest 
energy levels may have a contribution from the hot core 
(dense and hot inner regions). However, the large far-IR 
line-plus-continuum opacity will probably hide most of 
the hot core emission. Finally, radiative transfer effects 
in the most opaque lines, e.g. self-absorption and/or 
scattering by a lower density diffuse halo may possibly 
occur at velocity scales not resolved by ISO. 

In order to estimate the water abundance, and the 
physical conditions prevailing in the expanding gas, we 
have modeled the first 30 rotational levels of both ortho- 
and para-H20 using the same nonlocal code for lines 
and dus t continuum as in the analysis of H2O toward 
Sgr B2 IjCernicharo et al.ll2006tl. The code has been de- 
scribed elsewhere l|Gonzalez- Alfonso fc Cernicharoll993|) 
and has been recently improved to take into account a 
more sophisticated description of the dust emissivity and 
radiative transfer (Daniel et al. 2006, in prep.). The 
dust continuum emission has a crucial role in the ra- 
diative excitation of light species such as H2O or OH 
with transitions in the far-IR and has to be correctly 
taken into account. In our model, level populations 
are consistently computed in statistical equilibrium con- 
sidering collisional excitation and pumping by line and 
continuum photons. Collisi onal rates were s caled from 
those of H2 0-He collisions ijGreen et a]| ll993'). A three- 
component model resembling the hot core, the plateau 
and the ridge is adopted. In this work we assume that 
most of the far-IR water lines arise only from the plateau. 
Of course, a minor contribution in the lowest excitation 
water lines may come from the more extended regions 
(producing narrow line emission) that we don't model 
here. A low water abundance of (1-8) xlO -8 has been 
estimate d in these regions fr om SWAS and ODIN obser- 
vations l|Snell et alJ (|2000|k lOlofsson etafl l|2003|) 1 and 
thus, the expected contribution to far-IR H2O lines is 
small. In order to have the closest view of the IR radi- 
ation field seen by water vapor we have simultaneously 
tried to reproduce the full SWS/LWS/ISO continuum 
emission between -40 and -497 fim. 

We find that the continuum level for A >80 /j,m is 
well reproduced by considering an inner 10" IR contin- 
uum source (the hot core) simulated by a grey-body, 
with a color temperature of 200-250K, which is opti- 
cally thick in the far-IR (r A = 10 x 150/A). The dust 
continuum emission from this hot core is attenuated 
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Fig. 2. — Radiative transfer model results discussed in the 
text. Selected synthetic line profiles (continuous curves) are shown 
over the observed ISO/LWS/FP detections (histograms). Model 
parameters are x(H2 6 0)=3x 10 - 5 , n(H2)=3.5xl0 5 cm -3 and 
T fe =80 K (blue) and x(H2 6 0)=2xlCT 5 , n(H 2 )=2.5x 10 5 cm" 3 
and Tfc=100 K (red). The ordinate corresponds to the continuum 
normalized flux and the abscissa to the velocity (in km s~ 1 ). Syn- 
thetic line profiles have been convolved with a gaussian beam of 
80" and with a spectral resolution characterized by a Lorentzian 
with a width of 33 km s _1 . 



by the surrounding components that we now describe. 
Taking into account the spatial extent of the 183 and 
325 GH z water emission observed at higher spatial res- 
olution l|Cernicharo et al.lll990l 11994 H999bD , we model 
the plateau as a ~40" diameter unresolved shell expand- 
ing at 25±5 km s , where water molecules and dust 
grains coexist (Td=T\, is assumed). The radiative effects 
caused by an additional, only dust, 5' diameter compo- 
nent (the ridge), with T d ~25 K and n(H 2 )~10 5 cm~ 3 , 
have been also included to fully match the SWS/LWS 
continuum level. 

Dust absorption coefficients are computed from tabu- 
lated optical grain properties from Draine & Lee (1984), 
both in the plateau and in the ridge. Due to the large 
optical depth of H3, 6 lines in the plateau, we have first 
modeled the H^O lines to obtain a more accurate esti- 
mate of the water abundance with less opacity constraint 
and faster level population convergence. A 16 0/ 18 
abundance ratio of 500 has been adopted. After a rea- 
sonable fit to the H^O lines, we iteratively search for 
the best physical conditions that enable to simultane- 
ously reproduce the H^O and H^O spectra. Fig. 
shows two representative models that qualitatively and 
quantitatively reproduce the majority of observed line 
profiles and intensities. The best models are found for 
a plateau temperature, density and water vapor abun- 
dance of T fc =80-100 K, n(H 2 )=2.5-3.5xl0 5 cm' 3 and 



x(H2 6 0)=2-3x 10~ 5 respectively. The physical condi- 
tions inferred from water lines i n the plateau are simila r 
to those obtained from OH lines ijGoicoechea et ail 2006) . 
which suggests that far-IR lines from both species trace 
the same expanding gas. The derived water vapor abun- 
dance, >2xl0 -5 , is obviously an averaged value over the 
large LWS/ISO beam and probably indicates that water 
can be locally more abundant, e.g. in the warmer inter- 
action surfaces where the expanding gas shocks the dense 
ambient material. Larger angular resolution is needed to 
resolve the relative continuum and H2O spatial distri- 
bution over this complex region. Nevertheless, the in- 
ferred gas temperatures, ^100 K, are significantly below 
the gas temperature (~300 K) required to activate the 
gas-phase neutral-neutral reactions converting most of 
the available oxygen into water (abundances larger than 
~10 -4 are then predicted). Taken into account the short 
dynamical timescale of the outflow, ~1000 years, water 
vapor could have been formed by these neutral-neutral 
reactions if the gas in the plateau was previously (or is 
locally) heated to >500 K, e.g. by a C-shock passage 
ijBergin et alJ ll998L However, if the plateau gas tem- 
perature has reached a maximum temperature of only 
~100 K, other formation mechanisms are required to ex- 
plain the observed water vapor abundances. In partic- 
ular, the temperature inferred from ISO observations is 
similar to the water-ice evaporation temperature, and 
therefore, in-situ evaporation of water-ice grain mantles 
formed in an evolutionary stage prior to the onset of the 
outflow(s) could now dominate the water vapor forma- 
tion in the plateau. Finally, the observed H2O could have 
been also formed in the innermost and warmest regions 
(e.g. the hot core) and had been swept up by the outflow. 

To conclude, despite the large number of detected far- 
IR excited water lines, neither high gas temperature 
nor high density conditions are required to populate the 
higher-energy water rotational levels. Radiative pump- 
ing due to the strong IR radiation field from the cloud is 
enough to populate these levels, at least in the average 
picture given by ISO observations. Future observations 
with Herschel will allow one to map many far-IR water 
resolved lines with improved angular resolution. It will 
then be possible to study in great detail the role of water 
in star forming regions. 
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